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Abstract 

We investigated the bubble collisions during the first order phase tran- 
sitions. Numerical results indicate that within the certain range of param- 
eters the collision of two bubbles leads to formation of separate relatively 
long-lived quasilumps - configurations filled with scalar field oscillating 
around the true vacuum state. Energy is perfectly localized, and density 
is slightly pulsating around its maximum. This process is accompanied 
by radiation of scalar waves. 

I. Dymnikova, L. Kozicl Institute of Mathematics, Informatics and Physics, 
University of Olsztyn, Zolnicrska 14, 10-561 Olsztyn, Poland 

and M. Khlopov, S. Rubin Centre for Cosmoparticle Physics "Cosmion" 
Moscow State Engineering Physics Institute, Moscow, Russia 

PACS 04.70.Bw; 04.20.Dw 



1 Introduction 

Analysis of physical processes predicted in the early Universe on the basis of 
particle theory, is the important way to study physical conditions in the early 
Universe and physical mechanisms underlying those conditions. As a result of 
such an analysis, the existence of hypothetical relics of early Universe, such as 
primordial black holes, topologically stable or metastable solitons etc, have been 
predicted. Confrontation of predicted effects with observational data provides 
certain conclusions concerning both cosmological evolution and particle physics 
models 0. 

First order phase transitions as predicted by unified theories can occur at 
several periods of cosmological evolution. They are predicted on the basis of a 
wide class of models of particle symmetry breaking Q . They are also considered 
as the final stage of inflation in a wide class of inflationary models. Detailed 
study of nonlinear configurations arising at the first order phase transitions and 
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their dynamics is helpful not only for cosmology - nonlinear dynamics of field 
theories describes a wide range of phenomena occuring in laboratory physics. 

The transition with symmetry breaking consists in decay of a metastable 
phase by nucleation of bubbles of new phase |3| . A nucleated bubble is a true 
vacuum fluctuation large enough to evolve classically. The most likely fluctua- 
tion is a spherical bubble nucleated at rest with a certain critical size determined 
by microphysical processes j|. Coleman || calculated the bubble nucleation 
rate in flat space and at zero temperature using the euclidean path-integral for- 
mulation of a scalar field theory. The nucleation rate in this case is proportional 
to e~ E where Se is the euclidean action and solution to the euclidean equation 
of motion for minimal action is the 0(4) symmetric "bounce" solution. 

In the very early Universe phase transitions most likely occur at a finite 
temperature due to the fact that the form of a scalar field potential becomes 
temperature-dependent when quantum corrections are taken into account S. 
Generalization of Coleman results to the case of nonzero temperature is based on 
the remarkable fact that quantum statistics at nonzero temperature is formally 
equivalent to quantum field theory in the euclidean space, which is periodic 
in time coordinate with the period T -1 . As a result, most likely fluctuations 
appear to be not 0(4) symmetric spherical bubbles but 0(3) symmetric (with 
respect to spatial coordinates) cylindric configurations with certain critical size 
slightly different from 0(4) symmetric case [Q, 9. 

For bubble created with a size smaller than the critical one, it could seem 
that the gain in volume energy cannot compensate for the loss in surface en- 
ergy and such the bubbles would have to quickly shrink to nothing. However, 
detailed analysis discovered that even in this case effects of nonlinearity lead to 
nontrivial dynamics. The evolution of subcritical bubbles - unstable spherically 
symmetric solutions of nonlinear Klein-Gordon equation - was, firstly, studied 
numerically by Bogolubsky and Makhankov [|| . Using a quasiplanar initial con- 
figuration for the bubbles, they found that for a certain range of initial radii, 
the bubble, after radiating most of its initial energy, settled into long-lived (as 
compared with characteristic time-scale) stage and only then disappeared by 
quick radiating their remaining energy. Those configurations called "pulsons" 
were later rediscovered and revised by Gleiser who found that their most char- 
acteristic feature is not pulsating mechanism for radiating the initial energy, 
but the rapid oscillations of the amplitude of a scalar field during long-lived 
pseudo-stable regime when almost no energy was radiated away and radial pul- 
sations were rather small M . It was shown ^ that those configurations called 
"oscillons" exist for symmetric and asymmetric double- well potentials, are sta- 
ble against small radial perturbations, and have lifetimes "far exceeding naive 
expectations" ||. 

Although it is well known, that three-dimensional nontrivial configurations 
of a scalar field are unstable, they can be relevant for systems with short dy- 
namical time-scales. Detailed study of unstable but long-lived configurations 
can clarify dynamics of nonlinearities in field theories and their role in a wide 
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class of phenomena ranging from nonlinear optics to phase transitions both in 
the Universe and in the laboratory 0. 

For the bubbles formed with the radii large enough (overcritical bubbles) 
it is classically energetically favorable to grow. The newly formed bubble of 
true vacuum is separated from the surrounding false vacuum region by the wall 
at rest. Immediately after nucleation, the wall starts to accelerate outwards 
absorbing energy stored in false vacuum region and converting difference of 
false and true vacuum energy density into kinetic energy of the wall . That way 
a bubble spreads off converting false vacuum into the true one. This process 
continues up to the collision with a spherical wall of another bubble. 

In the first-order phase transitions at the end of inflation the collision of 
bubbles is considered as the leading mechanism of reheating by converting the 
wall energy into radiation. However situation with two bubbles appears much 
more complicated. Even nucleation of two bubbles is not yet studied in the 
literature in general [^) . Only in the case when bubbles are widely separated at 
the time of nucleation and thus can be treated as noninteracting (at the stage 
of nucleation) the generalization of a single bubble solution is straightforward. 

Two bubble collisions were studied in detail by Hawking, Moss, and Stewart 
[ flO| and then by Watkins and Widrow |J , in elegant approach using symmetry of 
the problem in zero temperature case. For zero temperature bubbles produced 
by quantum tunneling, initial state is 0(4) symmetric, as well as euclidean 
equation of motion, in natural assumption that a scalar field 4> is invariant under 
4-dimcnsional Euclidean rotations. In analytical continuation to Minkowski 
space this becomes 0(3, 1) symmetry. For two bubbles, the line joining their 
centers is the preferred axis and solution to the euclidean equation of motion is 
found in the class 0(3) (0(2, 1) as continued to the Minkowski space) solutions 
[ p"U[ and field configuration arising in collision belongs to the class of 0(3) 
symmetric solutions. 

The aim of present paper is investigation of two-bubble collision in the case 
of finite temperature. We are interested not in reheating by two bubble collisions 
[ pr] but in evolution of two bubble configuration during and after collision. The 
similar subject was discussed several years ago fll2fl . As was noted by Hawking, 
Moss and Stewart and confirmed by Watkins and Widrow, collision of two do- 
main walls does not lead to immediate convertion of the wall energy into a burst 
of radiation. Two walls reflect off one another and move apart creating a new re- 
gion of false vacuum between them |^] . Our aim was to investigate an evolution 
of this new false vacuum region to look if it can form a separated object. We 
connected with such a possibility the hope of formation of metastable relics of 
the first order phase transitions such as primordial black holes or selfgravitating 
particlelike structures with de Sitter- like cores G|. 

To our surprise, a false vacuum configuration evolved into a compact quasilump 
filled with an oscillating scalar field. The fundamental difference of this object 
from an oscillon is that it arises dynamically as the result of bubble collisions 
(which increases probability of its production) and that it is made up from an 
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oscillating scalar field at the background of true vacuum. We call it quasilump. 
since it does not satisfy all requirements for lumps as defined by S. Coleman: 
"Non-singular, non-dissipative solutions of finite energy, lumps of energy hold- 
ing themselves together by their own self-interactions" Jl4] |. In our numerical 
simulations we observe non-singular configurations of self-interacting scalar field 
with asymmetric potential, perfectly localized, but we cannot say that they are 
non-dissipative, although they are rather long-lived as compared with the char- 
acterstic scale for the first order phase transitions. 

Our paper is organized as follows. In Sect. 2 we present basic equations and 
initial configuration. In Sect. 3 we give qualitative analysis of the process of 
collision to find an optimal range of parameters for which concentration of false 
vacuum energy (which is later transformed into the energy of oscillations) is 
maximal. In Sect 4 we present the results of numerical simulations. Sect 5 
contains summary and discussion. 



2 Basic equations 

To study mechanism of formation and evolution of false vacuum regions, we 
shall consider most favorable regime for their appearance which corresponds to 
high nucleation rate NH A ^> 1 where N is the nucleation rate per unit four- 
volume and H is the Hubble parameter JlOj. We also neglect gravity effects 
on the process of bubble formation and growth which means that we consider 
bubbles with the initial size much less than cosmological horizon, R(0)H <C 1 
§§■ 

We consider real scalar field <j> with the Lagrangian 

c = ^d^cp-v(4>)- (i) 

This is the effective Lagrangian for a large number of more complex models of 
Universe involving the first order phase transitions (see Js) for more details). In 
the "thin wall" approximation, e/A << 1, some analytical results are known ||, 
and we will work in the frame of this approximation. 

To compare our results with the results obtained by Hawking, Moss and 
Stewart, and by Watkins and Widrow, we choose asymmetric double-well po- 
tential of the same form 

y(0) = iA(0 2 -^)2 + e0 3 (0 + 0o) _ (2) 

At T — the parameters A, <po and e are specified by the particle model. 
At nonzero temperature they are influenced by temperature corrections. In the 
case of high nucleation rate, a first order phase transition is a quick process 
and we can consider parameters as A ~A(T C ), 4>q ~ (j)o(T c ), e ~ e(T c ), i.e. being 
constant during the phase transition at the temperature T = T c . 
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The potential (|^) has two minima at different values of field (f> . False vacuum 
(metastable) state is characterized by the field <j> — ^o(l — e/A — 3/2(e/A) 2 ), 
whereas the global minimum of the potential V((f>) represents the true vacuum 
state cj) = -0 O (1 ~ e/A + 3/2(e/A) 2 ). 

In our analysis we assume that both mechanisms of the false vacuum decay 
could take place - tunneling, that is creation of 0(4) symmetrical bubbles, and 
formation of 0(3) symmetrical bubbles due to temperature fluctuation. Evi- 
dently, if the temperature is small enough, the tunneling mechanism of the false 
vacuum decay dominate. On the contrary, at large temperatures the decay is 
realized by the nucleation and growth of the 0(3) symmetrical bubbles. 

Consider conditions of dominance of the false vacuum decay due to temper- 
ature effects. The temperature decay probability was found in 0], ||: 

Premp oc e" S3/T , (3) 

where T is the temperature of a phase transition and S3 is three-dimensional 
action for 0(3) symmetrical bubble. The probability of the vacuum decay due 
to tunneling, is given by 

Ptun oc e~ Si (4) 

where S4 is the action for 0(4) symmetrical bubble. The temperature decay 
dominates if S3/T < S4. 

The straightforward calculations of the actions £3 and S4 give for our po- 
tential the condition for the dominance of the temperature decay (the term 
proportional to (e/A) 2 was omitted): 

in the units m v = 1. 

Consider the equation of motion of the scalar field in spherical coordinates 

d 2 cj> 9 2 2 64 (a , 
W ~ 9^ ~ = -V '(*), (6) 

Neglecting terms of order of 0((e/A) 2 ), we obtain the well known one- 
dimensional equation 

^-^-~ V W l-o • (7) 

The properties of this equation have been extensively discussed in the literature 
since 1975 flq]. The fundamental time independent solution is defined by 
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It can be easily checked by straight substitution, that for the theory defined by 
potential (||), the solution is represented in the form 

- MH^ir ~ R(f))] - e/X}, R(t) =vt + R , (8) 

where 7 = 1/^/1 — v 2 , v < 1, m = V\(f>o and Rq = 2A/(3em) is critical radius 
of the nucleated bubble. 

The initial field configuration can be defined at the moment of the bubble 
formation t = with velocity v = 0. But it takes too much computer time and 
memory to obtain the results of the collision, because the kinetic energy of the 
walls of the colliding bubbles should be large enough to produce false vacuum 
bag (FVB) and hence the initial distance between the centers of colliding bubbles 
should be large comparing with critical radius Rq as well. So, we have to use 
the initial configuration with already moving walls. 

The one-bubble solution (||) is the approximate solution to exact Eq.(^|). It 
also satisfies the correct boundary conditions at infinity up to the terms of order 
(e/A) 2 and hence can be chosen as new initial condition at definite moment t 
or at definite radius of the expanding bubble R = R(t). The only thing that 
remains to do is to connect the radius R and the velocity v. 

To find the velocity v in the one-bubble solution (||) at an arbitrary moment 
t or at definite bubble radius R(t) we note that the energy 

is conserved if the field 4> is governed by Eq. (^) . The substitution of the field <j) 
in form (|J) leads after simple calculations to the expression 

E~ ^ii?(t) 2 [ 7 -i?(i)e/A] = Const (9) 
3 A 

The Const can be determined at t = 0, because we know the values of 
the parameters at this moment: j(t = 0) = 1 and R(t = 0) = 2A/3e Q. 
Substituting it into expression (^|), we find the connection between the bubble 
radius R = R(t) and 7— factor (or, equivalently, the velocity v): 

^ R { + Ti^ ^ < 10 » 

Thus, the initial conditions for one bubble of radius R is represented by 
formula (H) with 7— factor (10). 



3 Bubble collisions - Qualitative analysis 

Let us introduce the dimensionlcss variables ip — </>/</>o, A 1 / 2 ^^ - * t and 
A 1//2 </>or — > r. The classical equation of motion for the scalar field of Lagrangian 
(Q) has the form 
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9 t V-V 2 V = -^W 2 -l)-e/A (11) 

The suitable initial two-bubble configuration has in our dimensionless variables 
the form 

^ = MHi{r+-R)]-e/\},z<Q, 
$ = M MUr- - R)} - e/A }, z > 0, (12 ) 
r± = \Jx 2 + y 2 + {z±bf, b > R. 

To start numerical simulations of bubble collisions, we need a set of sim- 
ple criteria indicating a proper range of parameters favorable to formation of 
separated false vacuum regions. 

Let us first find the condition at which the region of a false vacuum can be 
formed as a result of a collision of two relativistic bubbles. 

A field configuration in a wall is just transition from a true vacuum inside 
to a false vacuum outside. While propagating through a false vacuum before 
collision, bubble absorbs the energy of a surrounding false vacuum and trans- 
forms it into a kinetic energy of the wall. The kinetic energy is characterized by 
the Lorentz factor 7 = l/yl — v 2 . To get a region of a false vacuum between 
bubbles as a result of a collision, energy absorbed by walls from a false vacuum 
to the moment of a collision, must be sufficient to form a false vacuum state 
at least at the scale of the wall width. Let us estimate the lower limit for 7 at 
which such minimal region can be formed. 

Consider collision of two sphericall 0(3) bubble walls described by the so- 
lution (|^) with the parameters R and 7;„ to the moment of a collision. The 
leading term in the energy density of a wall, as calculated for the quasiplanar 
solution (||), is 

p w ~ 7 2 /4cosh 4 ( 7 (i?-^)/2) 
. Before a collision, in the solid angle 



A^-«l 

, each wall has the energy 

2 9 
E ln = -AOi? 2 7in . 

After a collision, the walls reflect with a final kinetic energy Efi n . If we want 
bubble wall collision to form, between reflecting walls, a false vacuum region of 
a radius r and width h within the solid angle AO, we must have 

2 2 £ 

E in — Efi n = -AOi? J f in + 2 — Vf vr , 
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where p vac = 2e/A is false vacuum density for the case of the potential (||), and 
Vf vr is the volume of a false vacuum region within a cone with a solid angle Afl 
which is given by 

Vf vr = ^ h (3r 2 + h 2 ) « ~AJXR 2 /i 

The width of a false vacuum region is of order of a width of a wall to the moment 
of reflection which is 2/jf in . For -ff in — 1, the width is h = 2. It gives us the 
constraint for ji n with which a wall came to the first collision in the form 

lin ~ 1 + ^J k - 1 + 3 I' ^ 

Now let us specify the line joining centers of bubbles as z axis. Let us show that 
the energy conservation puts constraint on the propagation of a false vacuum 
region in z direction. 

Consider a slice of a false vacuum region originated from the collision in 
the element Af2 of spherical bubble walls which has radius R in the moment 
of collision. The acceleration of the considered element of the wall came from 
transformation of the energy of surrounding false vacuum into kinetic energy of 
the wall on the way to its first collision, when the true vacuum bubbles grow 
from the initial radius R(0) to the radius R ^> R(0) in the moment of collision. 
So, the kinetic energy absorbed by the wall from a false vacuum to the moment 
of collision, is 

The walls reflect each other in the moment of the first collision and move out- 
wards, creating a false vacuum region between them. Each wall stops when all 
its kinetic energy has been transformed into the energy of a false vacuum region 
formed between the walls. In this moment the walls radius is R m ax and a false 
vacuum, created by each wall, fills a region between the spherical shells R max 
and R. The energy balance gives 

2e 1 o 2e 1 , , 
-AO-i? 3 ~ -Aft- (RL ax - R 3 ), 
A3 A 3 ' 

so that 

R max ~ 2 1 ^R. (14) 

Since we consider overcritical bubbles, the wall surface energy is neglected in 
this treatment, provided that ei? maa; /A ^> 1 The same result has been obtained 
for the case of 0(4) symmetric bubbles by HMS po| . 

One finds from the equation (|l4|) that after the collision a false vacuum is 
formed and occupates a region between the outgoing walls, with a maximal size 
given by a distance between the planes z = ±(2 1 / 3 — 1)6, where 2b is the initial 
separation of the centers of true vacuum bubbles. 
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After the walls stop their outward movement at R m ax — ±(2 1 / 3 — 1)R in 
the region of walls intersection, the parts of walls in this region reflect off one 
another and next time they collide at At ~ 2(2 1 / 3 — l)R after the first collision. 
The shortest interval between two subsequent collisions is at r — when R = b. 
Using the condition (|l^), we find from the Equation ( |lO| ) the minimal 7 at which 
a false vacuum region is formed between the walls after the second collision. 
Before the second collision at r = 0, the value of 7 in the walls is given by 

72 = R^k (2 1 / 3 - (15) 

Remember that before the first collision this factor for the walls at r = is 



given by the Eq.fllOD as 

71 = bj. (16) 

If we want the false vacuum region be maintained after the second collision of 
the reflected parts of walls, we must satisfy 72 > 7m- Then it follows from the 
Equations ( |13|JlqJlq ) , that 7 before the first collision must be 

l + 3e/A 

7l > 7mm - , 2 jy 3 _ ^ • (17) 

It indicates the favorable range of the 7 parameter before the first collision 
needed for numerical simulation, and also, by pC|Jl^ ), the favorable range for 
the parameters R and b. 

In the case 71 3> 7mm, a false vacuum region undergoes the succession of 
oscillations - expansions and contractions - along the z axis in the region confined 

by 

-6(2 1/3 + l) < z < 6(2 1/3 - 1). (18) 

Repeating the above reasoning for the subsequent collisions we find easily 
that in the limit of large 7 the period of the n— th oscillation decreases as 
(\f 7 2n)~ 1 . This agrees with the HMS result |l(J for the 0(4) symmetry case. 
The reason for such a coincidence can be easily understood. The main difference 
between the 0(3) and 0(4) cases is in the form of the initial wall configurations, 
taken in our case as a quasiplanar 0(3) solution. However, in all the above 
reasoning the internal structure of the walls was not involved which just resulted 
in similar estimation for the decreasing of period of oscillations. 

For large 7 we can treat the oscillations of a false vacuum region along the 
axis z as the continuous propagation of a spherical wave moving with speed of 
light (in our units c = 1). In the frame with the origin in, say, z — —6, the 
element Af2 of the wall with the angle a with respect to the z axis, follows 
the trajectory r = ztana. Assume that to the moment of reflection considered 
element of the wall has the coordinate z — ct. At the same moment its radial 
coordinate is r = ztana. The region of causal contact along the axis r satisfies 
the condition r < ct. It follows then that only for the angles a < 7r/4, the 
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region of intersection of walls is in causal contact. Therefore the boundary of 
the region of causal contact within a false vacuum region is the cone 



a = it/A. (19) 

It means that further evolution of a false vacuum region confined within the 
boundary ([l^) does occur independently on dynamics of field outside this bound- 
ary. So, considered region of a false vacuum is separated in its further causal 
and hence dynamical evolution. 

Now we can easily estimate the total energy of a separated false vacuum 
region. The energy density of a false vacuum is given by p vac = 2e/A. The 
volume of a FVB is the volume of two cones whose height is equal b and base 
area irb 2 . So, the mass confined within this region is 

4e irb 3 , . 

M = IT < 20 > 

It is evident that separation occurs at the time of order of t sep = (V2 — l)b. 



4 Numerical Results 

In the cylindric coordinates the equation of motion for a scalar field ( pd| ) equation 
takes the form 

dfy - a r V - 0ty = -\^ 2 - 1) - e/A 

The solution to this equation has the axial symmetry and reflection symmetry 
with respect to z = plane. The initial configuration described by solutions 

ip = f%/2(r + - R - vt)\ + tft[7/2(r_ -R- vt)] -1-e/X (21) 

is shown in Fig.l. The walls already have kinetic energy that is indicated by 7 
factor equals 5. 

Time evolution of a scalar field in the center of the region of collision ip(t,r = 
z = 0) shown in Fig. 2, was calculated for the parameters r = 5; b = 52; R = 50. 
The qualitative behavior of the field with time has been discussed in previous 
section. From the beginning the field changes in the manner discussed in JTof , 
then it oscillates around true vacuum for a long time and finally, large 
secondary fluctuations appear again. 

As we shall see below, the energy of oscillations is perfectly localized. This 
behaviour does not change with changing the step in calculations. Field config- 
urations in different moments of time are shown in Figs. 3, 4. 

The energy density profile calculated from scalar field potential, is shown in 
the next series of figures. They demonstrate concentration of field energy in the 
center of region of collision. In the fig. 5 one can see time dependence of energy 
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density in the center of the region of collision. Large secondary peak is created 
due to the coherent field oscillation which are coming from outside. 

The density profile at this time is represented in the Fig. 6. This localized 
configuration oscillates for some time and finally is converted into outgoing 
radiation. Only gravity could prevent this process. Till now we did not consider 
gravitational effects, but they will be estimated below. Consider the evolution 
of energy contained in the sphere of certain radius as shown in Fig. 7-8. 

This pictures show two peaks of energy - the first is due to energy in the 
moment of collision, the second is the energy of the quasilump which is formed 
as a result of the collision. It becomes evident by comparing these figures that 
the energy is strongly localized. Indeed, the energy contained in internal sphere 
of radius 5 is only in 1.5 times smaller than that contained in external sphere, 
while the ratio of their volumes is 8. Gravitational radius of the quasilump in 
our units is 




If the first order phase transition happens at the end of inflation, the mass 
m of inflaton field is rather large and m v /m pl ~ 10~ 5 . Substituting this value 
into (p^ ) and the value of energy E « 1000 obtained from fig. 7 one can easily 
find the condition when gravitational radius is comparable with the size of 
the quasilump, which can be taken as r$ ss 5in our dimcnsionlcss units. This 
condition is satisfied if coupling constant A ~ 10 -8 . For A < 10 -8 gravitational 
forces become essential and the probability of black hole formation grows up to 
unity when A tends to 10~ 8 . 



5 Summary 

In this paper we give qualitative arguments supported by numerical simulation 
for the existence of long-lived fluctuation that arise as a result of a collision of 
two expanded bubbles. 

The two-bubble collision leads, first, to the formation of short-living false 
vacuum region in the center of collision. Numerical results indicate separation 
of a false vacuum region at the time t ~ b. Then it evolves into rather compact 
object - quasilump made up of a scalar field oscillating around its true minimum, 
with lifetime enough to be captured by its gravitational field. At small coupling 
constants black hole can be produced. 

Till now the similar object discussed in literature was oscillon The 
main difference between these two objects is as follows, i) Oscillon represent 
a subcritical bubble of true vacuum inside a false vacuum, that arise due to 
temperature fluctuations. Our object is the fluctuation of scalar field in the 
true vacuum background that arise as a result of dynamical process, ii) To 
be long-lived, oscillon should have rather large initial radius, though less than 
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critical one, and rather flat initial distribution of scalar field. The evolution of 
the oscillon consists of oscillation of field value with almost constant radius of the 
field configuration. Our lump is much more compact object with the amplitude 
value of scalar field being much larger than that of the field in its potential 
minimum. For A < 10~ 8 gravitational forces are essential and the probability 
of PBH creation is of order unity, iii) Oscillon, being produced in spite of small 
probability fTi| , is extremely long lived object with lifetime 103 -104 1/m, m 
being the mass of scalar field. The life-time of our lump is of the same order of 
magnitude but lumps can be produced with much heigher probability, because 
they result from collisions of overcritical bubbles whose rate of nucleation is 
much bigger than for subcritical bubbles. 
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